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ABSTRACT 

Context. Low-mass star-forming cores differ from their surrounding molecular cloud in turbulence, shape, and density structure. 
Aims. We aim to understand how dense cores form out of the less dense cloud material by studying the connection between these two 
regimes. 

Methods. We observed the L1517 dark cloud in C 18 O(l-0), N 2 H + (J=l-0), and SO(J N =3 2 -2i) with the FCRAO 14m telescope, and 
in the 1.2mm dust continuum with the IRAM 30m telescope. 

Results. Most of the gas in the cloud lies in four filaments that have typical lengths of 0.5 pc. Five starless cores are embedded in these 
filaments and have chemical compositions indicative of different evolutionary stages. The filaments have radial profiles of C 18 0(1- 
0) emission with a central flattened region and a power-law tail, and can be fitted approximately as isothermal, pressure-supported 
cylinders. The filaments, in addition, are extremely quiescent. They have subsonic internal motions and are coherent in velocity over 
their whole length. The large-scale motions in the filaments can be used to predict the velocity inside the cores, indicating that core 
formation has not decoupled the dense gas kinematically from its parental material. In two filaments, these large-scale motions consist 
of oscillations in the velocity centroid, and a simple kinematic model suggests that they may be related to core-forming flows. 
Conclusions. Core formation in L1517 seems to have occurred in two steps. First, the subsonic, velocity-coherent filaments have 
condensed out of the more turbulent ambient cloud. Then, the cores fragmented quasi-statically and inherited the kinematics of the 
filaments. Turbulence dissipation has therefore occurred mostly on scales on the order of 0.5 pc or larger, and seems to have played a 
small role in the formation of the individual cores. 
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Dense cores in nearby dark clouds are the birth sites of solar- 
mass stars, and they represent the simplest environments where 
stars can be formed. The study of the internal structure of 
these cores offers a unique opportunity to determine the ini- 
tial conditions of low-mass star formation, and for this rea- 
son an intense observational effort is under way to characterize 
cores in detail (see recent reviews by di Francesco et alj|2007t 
IWard-Thomnson et alJl2007t iBergin & Tafallall2007l) . 

Compared to their turbulent parent clouds dense cores are 
very quiescent. Molecular line observations in tracers like am- 
monia reveal that the gas in a dense core has subsonic internal 
motions, indicating that thermal pressure contributes mo re tha n 
turbulence in providing support against gravity dMversl fl983). 
In addition, the gas motions inside a core are "coherent," in 
the sense that the ob served linewidth does not depend on radius 
(Goodma n et al.l 19981) . and thus deviates from the linewidth-size 
relati on characteristic of the large-scale gas in a cloud (Larson 
1 1981b . Also in contrast to the irregularly shaped cloud gas, cores 
often display a regular, close-to-spherical geometry ( Myers et al. 
1991). A number of cores even present centrally concentrated 
density profiles close to those predicted by hydrostatic equilib- 
rium models (Alves et al. 2001), again suggesting that the core 



gas conditions differ significantly from those of the cloud gas on 
larger scales. 

How quiescent, centrally-concentrated cores form out of 
the more turbulent, less-dense cloud material is still a mat- 
ter of debate. A number of core-formation mechanisms have 
been proposed over the years, ranging from the q uasi-static 
loss o f magnetic support via ambipolar diffusion (Sh u et all 
119871: iMouschovias & Ciolekl [T999) to the r apid dissipa- 
tion of turbulence due to supersonic shocks (Padoan et al. 
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2001; iKlessen et alJ l2005t IVazquez-Semadeni et all 12005). 
Observationally, it is unclear which of these proposed scenarios 
can fit the variety of existing constraints. Ambipolar diffusion 
models, for example, seem to predict stronger magnetic fields 
than observed (Crutcheret al. 2010), and they usually predict 
contraction times that significantly exceed thos e inferred from 
obse rved core life times and chemical clocks dLee et alJ ll999: 
Tafal la et al.ll2002l) . Models of core formation by shocks, on the 
other hand, predict linewidths and velocity displacements be- 
tween tracers that are larger than co mmonly observed in cores 
(IWalsh et alJl2004l: iKirk et alll2007a1) . 

Observational progress in constraining the mechanism of 
core formation requires probing the connection between the 
dense cores and the less-dense material that surrounds them. 
This less-dense material likely represents the gas out of which 
the cores have condensed, so by comparing its geometry, den- 
sity, and kinematics with those of the core gas, we could infer 
the physical changes involved in the formation of a dense core. 
Unfortunately, studying the cloud-to-core transition requires 
combining observations of tracers sensitive to different density 
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regimes, and until recently, such multi-tracer analysis has been 
hindered by a number of inconsistencies between the emission 
from different molecules (Zhou et al. 1989; Lemme et al. 1996). 

Work carried out over the past decade has revealed that most 
inconsistencies between tracers result from chemical changes 
that occur in the gas at dense-core densities and, in par- 
ticular, from the selective freeze-out of molecules onto cold 
dust grains (Kuiperet al. 1996; Caselli et al. 1999; Bergin et al. 
l200ilTafalla et al.ll2002l;fAikawa et alJliool . As a result of this 
work, a reasonably consistent picture of the chemical behavior 
of the different gas tracers has emerged. According to this pic- 
ture, most carbon-bearing molecules (including CO and CS) dis- 
appear from the gas phase at densities close to a few 10 4 cm" 3 , 
while nitrogen-bearing molecules like NH3 and N2H + survive 
undepleted up to densities that are at least one order of magni- 
tude higher. In addition, several species like SO and C2S present 
significant abundance enhancements at early evolutionary times, 
but end up freezing out onto the grains at simil ar densities to 
the C-bearing species (see Bergin & Tafalla 2007 for a review of 
core chemistry). 

This new understanding of molecular chemistry at interme- 
diate gas densities now makes it possible to combine observa- 
tions of different tracers and to finally explore the transition 
from cloud gas to core material self-consistently. An excellent 
region for investigating this transition is the L 1 5 1 7 dark cloud in 
the Taurus-Auriga molecular complex (see Ke n von et al.1 2008 
for a Taurus overview). L1517 appears in optical images as 
a region of enhanced obscuration associated with the reflec- 
tion nebulos i ty from the young stars AB Aur and SU Au r 
(iLvndsll 19621: IStrometal.lll976t [Schneider & Elmegreenlll979l) . 
The CO observations of the cloud have shown that the extended 
gas consists of several filamentary components that extend to 
the northwest of the nebulosity and occupy a regio n of about 
20' x 10' coincident with the optical obscuration (Heveret al. 
1987). Embedded in these components lies a collection of sev- 
eral dense cores that are easily distinguished in the optical im- 
ages thanks to the contrast provided by the bright nebulosity 
from AB Aur and SU Aur, which lie at least 0.3 pc in projection 
from the cores ( Schneider & Elmegreen 1979). Radio observa- 
tions of these cores reveal the phy sical properties typica l of the 
Taurus- Auriga core population (IBenson & Mverslll989l) . and a 
prominent member of the group, L1517B, has been the subject 
of detailed study owing to its regular shape and clear pattern of 
molecular freeze out dTafalla et al. 2004. 20061). 

A notable featu re of the L15 17 cores is that they a l l appear to 
be star less dStrometalJll976t lBeichmanet~aT1ll986t iKirk et alJ 
l2007bl) . This lack of embedded protostars is probably responsi- 
ble for the L1517 cloud presenting some of the narrowest line 
profiles in CO and dense gas tracers seen towards the Taurus- 
Auriga region, and indicates that, although the bright PMS star 
AB Aur (~ 50 L Q , van den Anck eret al][l9~9 8) is physically re- 
lated to the cloud (Nachman 1979; Du vert et aD fl986). the bulk 
of the L1517 material and its embedded cores remain unper- 
turbed by the energetic stellar output ( L add & Mversll99ll) . This 
combination of quiescent state, compact size, and multiple star- 
less cores makes the L1517 cloud an ideal laboratory for study- 
ing the process of core formation. In this paper, we present a 
study of the physical conditions and kinematics in both the dense 
cores and the less-dense surrounding material using a variety of 
molecules known to trace different density regimes and chem- 
ical evolutionary stages of the cloud gas. As will be seen, our 
analysis suggests that the cores in L1517 have formed by the 
gravitational contraction of subsonic, velocity coherent gas in 
0.5 pc-long filaments. 



2. Observations 

We observed the L1517 dark cloud with the FCRAO 14m ra- 
dio telescope during several sessions between December 2003 
and November 2005. We used the 32-pixel SEQUOIA array to 
cover the cloud with five sub maps, each of them of 10' x 10', and 
observed in on-the-fly mode. The large bandpass of SEQUOIA 
allowed observing two different transitions simultaneously, and 
two passes were made to cover the cloud first in N2H + (J=l-0) 
and SO(J N =3 2 -2i) and then in C 18 O(J=l-0) and C 17 O(J=l-0). 
All observations were made in position-switching mode using 
a reference position offset by (18', -13') from our map cen- 
ter (at a(/2000) = 4 A 55 m 18 s .8, 5(72000) = +30°38'04"), and 
known to have negligible C 18 O(J=l-0) emission from previ- 
ous frequency-switched observations (average cr(T m b) < 0.1 K 
over the SEQUOIA footprint). The spectrometer was the DCC 
autocorrelator configured to provide 1024 spectral channels of 
25 kHz spacing, or approximately 0.07 km s at the observ- 
ing frequencies. The telescope FWHM beam size varied with 
frequency between approximately 56" at the lowest (N2FU) fre- 
quency and 47" at the highest (C lv O) frequency. 

During the observations, calibration was achieved by mea- 
suring the emission from the sky and an ambient load every 
a; 10 minutes, and the derived intensity was converted into the 
main beam brightness temperature scale using facility-provided 
main beam efficiencies close to 0.5. The telescope pointing was 
checked and corrected approximately every three hours by mak- 
ing five-point maps on the SiO masers of IK Tau and Orion-IRc2. 
Typical pointing errors were within 5" rms. 

The off-line data reduction consisted in the creation of 
Nyquist-sampled maps with the otf tool program and a conver- 
sion to the CLASS format for further analysis with the GILDAS 
software (http : //www . iram. fr/IRAMFR/GILDAS). This anal- 
ysis included a second-degree baseline subtraction and a spatial 
convolution with a Gaussian to eliminate residual noise. The fi- 
nal resolution of all the FCRAO data is 60", taking the off-line 
convolution into account. 

Because the on-the-fly FCRAO maps are relatively shal- 
low because of the large area covered, we complemented the 
N2H + (l-0) FCRAO observations with a small amount of data 
observed with the IRAM 30m radio telescope. These data were 
taken in frequency-switching mode and have a higher angular 
resolution than the FCRAO data (about 26"), but apart from that, 
they are consistent with the FCRAO observations in both inten- 
sity and frequency calibration. A full report of these IRAM 30m 
data will be presented elsewhere as part of an extended study 
of core chemistry that includes a number of regions in addition 
to L1517 (Tafalla et al. in preparation). In the present paper, the 
IRAM 30m data are only used to complement the FCRAO ob- 
servations in the velocity analysis, which requires a high S/N 
threshold. 

Given the narrow lines and small velocity variations mea- 
sured towards L1517, accurate rest-line frequencies were re- 
quired for analyzing of the spectra. For this reason, we 
used the most recent laboratory estimates for C 18 O(J=l-0) 
(109782.176 MHz, see lCazzoli et a l. 2003) and for C lv O(J=l- 
0) (1123 58.990 MHz for the brightest, J F=l 7/2-0 5/2 compo- 
nent, see lCazzoli et al.ll2.002h . which have 1 kHz or better accu- 
racy. For N 2 H + and SO, no accurate laboratory measurements 
are av ailable, so we used astronomical estimates. Following 
Pagani et al. (2009), we assumed a frequency of 93173.764MHz 
for N2H + (Jf iF = 123 -012), with a n estim ated uncertainty of 
4 kHz, and following Tafall aetal] (|2P06), we assumed a fre- 
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Fig. 1. Map of the C 18 O(l-0) integrated intensity towards the L1517 cloud superposed to the red DSS image. The emission has 
been integrated over the whole velocity range of detection (Vlsr = 5.25-6.45 km s -1 ). First contour and contour interval are 
0.22 K km s . The two bright stars surrounded by nebulosity and located near the southeastern edge of the C 18 O(l-0) emission are 
AB Aur (brightest) and SU Aur. 



quency of 99299.890 MHz for SO(J N =3 2 -2i), with an estimated 
uncertainty of 10 kHz. 

Additional observations of the dense cores in L1517 were 
carried out in the 1 .2mm continuum using the MAMBO array on 
the IRAM 30m radio telescope during several sessions between 
December 1999 and January 2 005. The data for co re L1517B 
have already been presented in Tafall a etalj (l2004h . while the 
data for cores L1517A and L1517C are newly reported here. 
(L1517D was not observed.) In all cases, the observations were 
done in on-the-fly mode scanning the telescope in azimuth at a 
speed of 4"s _1 , and using wobbler switching with a period of 
0.5s and a throw of 53" or 70". The raw data were corrected 
for atmospheric attenuation using sky dips usually made before 
and/or after the source observation, and the absolute intensity 
calibration was achieved using facility-provided factors derived 
from observations of planets. The older L1517B data were re- 
duced with the NIC software using no method of noise reduc- 
tion (to avoid filtering the extended emission), while the newer 
L1517A and L1517C data were reduced with MOPSIC using a 
noise reduction method optimized to recover the extended emis- 
sion. The intrinsic beam size of all the continuum data is 11", 
although the maps were later convolved to an equivalent resolu- 
tion of 20" to eliminate high spatial-frequency noise. 



3. CO data: the extended cloud 

Figure Q] presents our map of C I8 O(l-0) emission integrated 
over the full velocity range of detection superposed to the red 
DSS ima ge. This map is in good agreement with the previous 
map by iHever et al] dl987l) . and it shows how the C 18 O(l-0) 
emission extends to the northwest of AB and SU Aur, the bright- 
est stars in the DSS image. These PMS stars are physically as- 
sociated with the LI 5 17 cloud, as indicated by the heatin g of 
the gas in their vicinity (Nachman 1979c iDuvert et al.l il986). but 



their effect must be highly localized since the data described in 
the next sections show no evidence of any kinematic interaction 
between the gas and the stars. This lack of kinematical inter- 
action, together with the fact that the nort hwest elongation of 
the L1517 cloud continues in the maps of Duv ert et al.l (1986) 
to scales as large as 4pc (or more than 5 times the size of our 
map) and involves additional dark clouds such as L1496, L1505, 
L1513, and L1515, suggests that the distribution of gas in L1517 
is intrinsic to the cloud, and that it has not been sculpted by the 
PMS stars . In this sense, the L 1 5 1 7 cloud seems a quiescent rem- 
nant of a larger gas cloud that gave rise to AB Aur, SU Aur, 
and an additional gr oup of PMS stars located in their vicinity 
(Luhm anetail20 09). 



3.1. Filament identification 

A more detailed view of the gas distribution in L1517 comes 
from the analysis of its velocity structure. Figure [2] shows maps 
of C 18 O(l-0) emission integrated every 0.4 km s -1 , or about 
twice the sound speed for gas at 10K. These maps, and others 
made using different velocity ranges, show that the L1517 cloud 
is structured in at least four elongated components that are re- 
ferred to as filaments 1 to 4. These four filaments can also be 
seen in maps o f extinction (Jouni Kainulainen, private communi- 
cation), 13 CO (lHeveretal.ll 19871) . and SO emission (Sect.®, so 
they must reflect the true distribution of gas in the L1517 cloud, 
and are not mere artifacts of the C ls O chemistry or excitation. 

As Fig. [2] shows, each filament appears in only one 
0.4 km s~'-wide channel map, indicating that its gas is highly 
confined in velocity space. This low level of velocity structure 
results from a very quiescent state of the gas in the filaments, 
both in terms of velocity dispersion and internal velocity gra- 
dients (to be discussed below), and makes it possible to asso- 
ciate each element of emission in the cloud to one of the four 
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Fig. 2. Maps of C 18 O(l-0) intensity integrated every 0.4 km s ' to illustrate the separation of the gas into filaments. The dashed 
lines indicate our best guess of each filament axis, and the star symbols mark the positions of AB Aur (larger) and SU Aur (smaller). 
First contour and interval are 0.2 K km s _1 . Offsets are referred to a{J20QQ) = 4*55 m 18 8 .8, 5{J20QQ) = +30°38'04". The LSR 
velocity range of integration is indicated in the upper left corner of each panel. 



filaments. To carry out this assignment, we have taken both the 
spatial location of the emission and its velocity into account and 
assumed that the velocity limits of the filaments are those used 
in the maps of Fig. [2] (see also Table [TJ. This assigning proce- 
dure works well over most of the mapped region, but it provides 
ambiguous answers in the vicinity of offsets (-100", -300"). In 
this region, filaments 1 and 2 overlap spatially and converge in 
velocity to a value close to 5.65 km s , which is our assumed 
boundary between the two filaments. A detailed inspection of 
the FCRAO spectra in this region (supplemented with additional 
data from the IRAM 30m telescope) reveals that the emission 
from the two filaments becomes almost indistinguishable, either 
because the filaments merge physically or because they super- 
pose along the line of sight with similar kinematic properties. 
As a result, assigning emission in this region to either filament 1 
or 2 becomes more uncertain than in other places of the cloud. 
Although this may lead to some confusion in a localized region, 
it is unlikely to affect our global analysis of the gas kinematics, 
since the properties of the two filaments become so similar in 
terms of velocity centroid and linewidth that an erroneous as- 
signment of a gas parcel to either filament will only add a small 
amount of mass, but will not change its velocity field. 

Once we have decomposed the L1517 cloud into its con- 
stituent filaments, we use the C ls O (1-0) emission to estimate 
their basic physical parameters, and we summarize the results in 
Table Q] From the maps of Fig [2] we estimate that the filaments 
have lengths of 550-1000 arcsec, which correspond to physical 
sizes of ~ 0.35-0.70 pc for our adopted dista nce of 144 pc (based 
on th e Hipparcos distance to AB Aur, Ivan den Ancker et alj 
1998). These sizes represent a significant fraction of the total 
cloud length, and this reinforces the idea that the cloud is struc- 
tured as a network of filaments. To estimate the masses of the 
filaments, we use the C 18 (1-0) emission and assume that it 
is optically thi n and in LTE at 10 K. For a C 18 abundance 
of 1.7 x 10~ 7 dFrerkinp et alj |1982). the typical filament masses 
are 5-11 M G , which imply linear mass densities in the range 12- 
HMopc- 1 . 



3.2. Density structure of the filaments 

As Fig.|2]shows, the C ls O emission from the filaments displays 
a significant degree of central brightening. This suggests that the 
underlying density structure of the filaments is centrally con- 
centrated and that the C 18 emission can be used to determine 



Table 1. Filaments in L1517. 



Filament 



Vlsr range 
(kms- 1 ) 



Length'" 
(pc) 



Mass (2) 
(M ) 



Cores 



1 


[5.25 - 5.65] 


0.52 


8.0 


A2, C 


2 


[5.65 - 6.05] 


0.42 


7.2 


A1,B 


3 


[5.65 - 6.05] 


0.70 


11.3 


D 


4 


[6.05 - 6.45] 


0.38 


4.8 


— 



Notes: (1) Uncorrected for projection effects; (2) From C 18 O(l-0) 
emission. 



each filament density profile. To carry out this determination, we 
first defined a central axis for each filament assuming that they 
are rectilinear and show the results in Fig. [2] As can be seen, the 
rectilinear approximation is reasonable for all filaments but num- 
ber 3, which bends towards the location of the PMS stars near 
its southeast end. Whether this bend is related to core L1517D, 
with which it coincides in position (Sect. 4), or is associated to 
the PMS stars is unclear from our data. In either case, for sim- 
plicity in the modeling, we approximate all filaments by straight 
lines, with the caveat that this approach is only a first-order ap- 
proximation. 

Using the above linear axes, we created a profile of C I8 (1- 
0) intensity for each filament as a function of cylindrical radius. 
The result, presented in Fig. [3] shows that the four filaments in 
the cloud have centrally concentrated distributions with a flat- 
tened region close to the axis and an approximate power-law tail 
at large distances. Filament 3 presents the highest dispersion of 
the sample due to a combination of contamination from filament 
2 emission near its northern edge and the additional contribu- 
tion from the "bend" region previously discussed. Nevertheless, 
given the rather irregular distribution of emission seen in Fig. [2] 
it is remarkable that the radial profile of each filament presents 
relatively little scatter and that all the filaments seem to follow a 
very similar type of radial profile. 

The regular behavior of the radial profiles in Fig. [3] sug- 
gests that the underlying density structure of the filaments can 
be described with a simple density law. To derive such a law, 
we model the C ls O (1-0) emission assuming that each fil- 
ament has cylindrical symmetry and that the C 18 emission 
arises from gas in LTE thermalized at 10 K, as suggested 
by the study of the L1517B core and its surrounding enve- 
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Fig. 3. C 18 O(l-0) emission profiles across each of the L15 17 fil- 
aments. The points represent the integrated emission within the 
velocity range of each filament (see Table[T), and R/,7 is the dis- 
tance to the filament axis. For each filament, the red solid curve 
indicates the isothermal cylinder model that fits the emission to- 
wards the axis, and the blue dashed line is the best-fit softened 
power law. 



lope in Tafall a etaD (12.004ft . Under these conditions, if we as- 
sume a density pr ofile, together with a constant C 18 abun- 
dance (1.7xl(T 7 , Frerki ng et alll 19821) and a C I8 linewidth of 
0.3 km s" 1 (Sect. 16. il l, we can predict a radial profile of emission 



that can be compared, after appropriate beam convolution, with 
the observed radial profiles shown in Fig. [3] 

As a first guess for the filament density profiles, we use 
the family of isothermal cylinders in pressure equilibrium with 
their self gravity. T his family w as originally described by 
IStodolkiewiezl (1 19631) and lOstrikeri (fl964). and has the simple 
analytic form 



n(r) = 



«0 



(1 + (r/H) 7 -) 



2^2' 



where r is the cylindrical radius, n is the central density, 



H 2 = 



2d 



nGnn ' 



(1) 



(2) 



c s is the isothermal sound speed, G the gravitational constant, 
and fi the mean molecular mass. For the case of L1517, the 
gas kinetic temperature is known to be approximately 10 K 
(Tafalla et al. 2004) so any isothermal cylinder is described by 
just one parameter, its central density n . This means that once 
we have set the central density to a value that fits the intensity 
at zero radius, we do not have any free parameter left to control 
the rest of the radial profile, and in particular, the width at half 
maximum of the emission is automatically determined by the H 
value. 

There are a number of motivations for using the isothermal 
cylinder family as a first choice to fit the C 18 radial profiles. In 
addition to the inherent simplicity and assumption of equilibrium 
of the model, the isothermal cylinder is attractive for presenting 
both a central density flattening and a power-law tail, which are 
two characteristics of the observed radial profiles. The isother- 
mal cylinder, in addition, has a mass per unit length that is inde- 
pendent of t he central dens ity and is only a function of the gas 
temperature (Ostrik erl 19641) . For the 10 K assumed for the gas in 
L1517, this mass per unit length is approximately 16 M pc _1 , 
which is reasonably close to the mass per unit length of the fila- 
ments estimated in Table [1] 

The results of our isothermal cylinder fits to the filament pro- 
files in L1517 are shown in Fig. [3] These fits were selected to 
match the C 18 emission towards the filament axis, so the radius 
of half maximum emission is a direct prediction from the model 
so is not controlled by our fit. As can be seen, the predicted fila- 
ment widths are in reasonable agreement with the observations, 
with only filament 4 being clearly narrower than predicted by 
the model. In addition, the central densities predicted by the fits 
range from 3xl0 3 to lxlO 4 cirT 3 (see Tabled, and these values 
agree with previous esti mates of the gas den sity in the extended 
part of the L15 17 cloud dTafalla et alfeooH) . 

To check the consistency of the parameters derived with the 
isothermal cylinder model and to provide a better match to the 
radial profiles of the filaments, we have fitted the C ls O emission 
with an alternative family of profiles that also have central flat- 
tening and asymptotic power-law behavior. This family consists 
of softened power laws described by 



n(r) 



l + (r/r V2 ) a 



(3) 



where n is again the central density, ri/2 the half-density radius, 
and a the asymptotic power index. This family of profiles has 
been shown to fit the density structure of starless cores, which 
also present flatten ed central re gions and asymptotic power- 
law behavior (Tafal laet al.ll2004l) . In contrast to the equilibrium 
cylinder, the softened power law has three free parameters, so 
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it fits the central density and the width of the filament indepen- 
dently. By comparing the results from the softened power law fits 
with those from the isothermal cylinders, we can now test how 
self consistent the one-parameter isothermal-cylinder fits are. 

Figure [3] shows the softened power law fits derived with a 
chi-squared minimization algorithm. As can be seen, these fits 
are almost indistinguishable from the isothermal cylinder fits for 
filaments 1, 2, and 3, and the only clear difference between the 
two families of fits occurs in filament 4. A more quantitative 
comparison between the fits comes from examining the derived 
central density and filament width in the two families. Table [2] 
shows that these parameters differ by less than 10 % on average 
for the first three filaments, indicating again that the isothermal 
cylinder and softened power-law fits are equivalent within the 
scatter of the observations. In filament 4, the softened power- 
law fit is clearly superior to the (poor) isothermal cylinder fit, 
and the difference between the derived parameters in the density 
law is approximately a factor of 2. 



3.3. Implications and limitations of the modeling 

Our analysis of the L1517 filaments adds to a number of previ- 
ous stu dies of filamentary stru ctures. With few exceptions, like 
that of iJohnstone et alJ (120031) who fitted the 850 /mi emission 
from the infrared-dark cloud Gil. 11 -0.1 2, most previous work 
has found significant deviations between the radial profiles of 
the filaments and the prediction from the isothermal cylinder 
model, usually because the observed asymptotic power-law be- 
havior is flatter than the r~ 4 predicted by the isothermal cylinder 
mod el dAlvesetal.fi 998: Ste pnik et al.ll2003blAndre e7aTll2010t 
Arzoumanian et al. 2011). For the filaments in L1517, our ra- 
dial profiles do not extend far enough in radius to sample the r"" 4 
asymptotic behavior, as can be seen in Fig.[3]from the isothermal 
cylinder models being practically indistinguishable from soft- 
ened power-law profiles with an asymptotic behavior close to 
r~ 27 (Table|2]). As a result, testing whether the L1517 filaments 
follow the prediction from the model or deviate from it cannot be 
done with our data, and it requires extending the radial profiles 
to larger radii, preferably using more robust techniques, such as 
extinction measurements or dust continuum observations. 

Even if the L1517 filaments do not follow the expected 
asymptotic behavior, it is striking that they approximately fit the 
expectations from the isothermal cylinder model both in width 
and mass per unit length. Clearly there are a number of limi- 
tations in our study, like the use of the freeze-out prone C ls O 
molecule as a column density tracer and our ignoring of pro- 
jection effects, so there is probably room for a factor of 2 un- 
certainty in the results. Also, as we will see in section [8j the 
size scale of fragmentation suggests that there may be a con- 
tribution from additional forces, like external pressure or mag- 
netic fields (which cou ld even flatten the radial profile, e.g., 
iFiege & Pudritzl 1200 0a). Still, the very quiescent state of the 
cloud described in the following sections indicates that the gas 
in L1517 cannot be too far from a state of equilibrium, as 
otherwise it would quickly develop supersonic motions, (e.g., 
Burkert & Hartmann 2004) and would not fragment into well- 
separated cores (Inutsu ka & Mivamalll997l) . As pressure forces 
dominate, it seems reasonable to expect that the distribution of 
mass in the filaments bears some similarity to the prediction 
from the isothermal cylinder model. 



Table 2. Fits to the C 18 radial profiles 



(i) 



Isothermal cylinder 



Softened power law 



Filament 



(cm' 3 ) 



r < 2 ) 

(arcsec) 



(cm" 3 ) 



r 1/2 a 

(arcsec) 



1 


lxlO 4 


89 


8.5 x 10 3 


105 


3.4 


2 


7xl0 3 


106 


7xl0 3 


100 


2.7 


3 


8xl0 3 


99 


8.5 x 10 3 


90 


2.7 


4 


3xl0 3 


162 


6.5 x 10 3 


75 


2.8 



Notes: (1) Not corrected for inclination; (2) ri /2 = (2 



1/2 



l)" 2 H 



4. Dense core population 

The popul ation of dense cores in the L15 17 cloud was first de- 
scribed by ISchneider & E lmegreenl dl979h . who identified four 
regions of enhanced obscuration in the Palomar optical plates 
and named them A to D in order of increasing right ascension 
(RA). For unknown reasons, this original labeling scheme was 
later replaced by one in which the cores are named C, A, B, and 
D in order of increasing RA, and t his notation has been us ed by 
most core studies in the past (e.g. jBenson &Mverslll989l) . For 
consistency with previous work, we also use this disordered core 
labeling here. 

Figure [4] presents different views of the core population in 
L1517. The optical DSS image on the top left has had its con- 
trast enhanced to better show the cores as dark patches against 
both the galactic stellar background and the bright diffuse emis- 
sion from AB Aur and SU Aur, which are probably slightly 
behind the cores. The other panels show maps of the same re- 
gion in the 1.2mm dust continuum, N2H + (l-0), and SO(32-2i), 
three tracers that highlight different properties of the dense, high 
column density gas. The 1.2mm continuum emission is mostly 
sensitive to the column density of the cores (the extended emis- 
sion from the filament is filtered out by the bolometer), and its 
maps provide high angular resolution views of cores A, B, and C 
(core D was not mapped). The brightest 1.2mm emission corre- 
sponds to core B, whose structure and chem i cal c omposition has 
been studied in detail bv lTafalla et~aTI (12004 [2006h . Cores A and 
C have received less attention due to their w eaker emission in 
both high-density tracers and mm-continuum ( Benso n& Mversl 
1989; Ladd & Myers 1991; Kirk et al. 2005), so their structure is 
less well known. Our mm-continuum map shows that core A is 
double-peaked, while core C is single-peaked but relatively more 
extended and significantly weaker than cores A and B. Overall, 
the mm-continuum map is characterized by a lack of point-like 
components, which together with the absence of Spitzer 24 /mi 
point sources (Kir k et aDl2007bl) is a strong indication that all 
dense cores in the cloud are starless. 

The N2H + and SO maps in Fig. |4] present complemen- 
tary views of the L1517 cores. N2H + is a so-called late-time 
molecule, and i ts abundance is furthe r enhanced when CO 
freezes out (e.g. jBergin & Tafalla 2007), so bright N2H + emis- 
sion is commonly associated with evolved cores dCrapsi et alJ 
2005). As Fig. |4] shows, the N2H + emission is brightest towards 
cores B and D, and still noticeable towards core Al. These cores 
therefore seem more chemically evolved than cores C and A2, 
which are barely detected in N2H + (l-0). Such interpretation of 
the N2H + maps in terms of chemical evolution is supported by 
the distribution of SO emission. Observation ally, SO is known t o 
be highly sensitive to molecular freeze out (Tafalla et al. 2006), 
and in addition, it is predicted to decrease in abundance with 
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Fig. 4. Core population of L1517. From left to right and from top to bottom: contrast-enhanced DSS blue image identifying the 
core positions, 1.2mm dust continuum emission, N2H + (l-0) integrated intensity, and SO(32-2i) integrated intensity. Offsets and 
star symbols are as in Fig. [2] First contour and spacing are 2 mJy/ll"-beam for 1.2mm continuum, 0.16 K km s~' for N2H + (l-0), 
and 0.08 K km s~' for SO(32-2i). To enhance the sensitivity to extended emission, the 1.2mm continuum map has been convolved 
to an equivalent resolution of 30", and the N2H + and SO maps to a resolution of 75". 



Table 3. Kinematic properties of the L1517 cores. 



Core 


Aa, AS 


V,,„-(N 2 H + ) 


AV(N 2 H + ) 


<|VV, [r (N 2 H + )|> 




(". ") 


(kms- 1 ) 


(kms- 1 ) 


(km s _1 pc -1 ) 


Al 


-60, -300 


5.70 


0.18 


0.7+0.1 


A2« 


-150, -240 


5.57 


0.20 


(2) 


B 


-30, -30 


5.79 


0.22 


0.7±0.3 


C (D 


-480, -150 


5.48 


0.18 


(2) 


D 


540,0 


5.88 


0.29 


0.9±0.4 



Notes: "' Kinematics data from IRAM 30m observations (rest from 
FCRAO 14m). <2) Not enough data to estimate gradient. 



time dBergin & Langerll997blAikawa et alJJOOlh . The SO emis- 
sion, therefore, is expected to be anticorrelated with that of 
N2I-F, and indeed, the maps of L1517 show such behavior: the 
SO emission is dominated by cores C and A2, which are weak in 
N 2 H + (l-0), while the N 2 H + -bright cores Al, B, and D are barely 
distinguishable from the SO emission of the extended cloud. 



4. 1 . Core emission modeling 

To quantify the physical and chemical properties of the L1517 
cores, we modeled the i r emis sion following the procedure de- 
scribed in lTafallaetai] d2004l) for the analysis of the L1498 and 
L1517B cores. In this way, we assumed that the cores are spher- 
ically symmetric, and we concentrated our modeling effort on 
fitting the radial profiles of emission shown in Fig. [5] To deter- 
mine the density profile of each core, we fitted the 1.2 mm con- 
tinuum emission, as this is expected to be the most faithful tracer 
of the total dust and gas column densities (e.g. jBergin & T afalla 
2007). Following the analysis of L1517B, we assumed a uni- 
form dust temperature of Tj - 10 K and a 1.2mm emissiv- 
ity of k — 0.005 cm 2 g , and we fitted the continuum radial 
profiles with density laws of the form n{r) - «o/(l + (r/ro) a ), 
where no, ro, and a are free parameters. (For core B, we have 
used the Tafalla et al. 2004 result, while core D was not fitted 
for lack of continuum data.) The results of these fits are il- 
lustrated in Fig. [5] and the best-fit parameters are summarized 
in Table [4] As can be seen, core central densities range from 
4.7 x 10 4 cm -3 in core C to 2.2 x 10 5 cm -3 in core B, which cor- 
respond, respectively, to enhancements of 6 and 30 with respect 
to the central density of the filaments traced in C 18 (section 
13.21 1. Integrating the density profiles up to a representative radius 
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of 0.05 pc (x 75"), we estimate that the core masses are on the 
order of 1-2 M a , which is typical of the population of Taurus 
starless cores (Benson & Myers 1989). 

Once the density profile of each core had been modeled, 
we could determine the abundance of the different species by 
solving the equation of radiative transfer and convolving the 
result with the appropriate gaussian beam size to simulate an 
observation. We did this using a Monte Carlo non-LTE code 
based on that of lBernesl d!979l) . and we excluded core D from 
the analysis due to the lack of mm-continuum data. For all the 
cores, we assumed a gas kinetic tem perature of 9 .5 K based 
on the NH 3 analysis of L1517B in Tafa lla et all ([2004), and 
in agreeme nt with typical es timates for other low-mass star- 
less cores dBenson & Mversl Il989t) . To reproduce the typical 
N 2 H + (l-0) FWHM of 0.21 km s _1 (mean of values for all cores, 
see Table [3]), we added a constant nonthermal component of 
0.17 km s (FWHM) to the velocity field. Also, we extended 
each core density profile with a cloud component based on the 
profile of the filament in which it is embedded. The assignment 
of cores to filaments was done using the core LSR velocity de- 
termined from N2H + (Table [3} and the velocity limits of each 
filament measured from C 18 (Table [TJ. According to this cri- 
terion, cores A2 and C belong to filament 1, cores Al and B to 
filament 2, and core D to filament 3 (as summarized in Table[TJ. 
As discussed in section 13.11 there is some ambiguity in the as- 
signment of gas to filaments in the region where filaments 1 and 
2 overlap, and this may affect the filament assignment for cores 
Al and A2. Our choice denotes our best understanding of the 
emission after a careful inspection of the individual spectra from 
both the FCRAO and IRAM 30m telescopes. 

The above parametrization fixes the physical structure of 
each core and leaves the abundance law of each molecular 
species as the only free parameter left to fit the observed radial 
profiles of intensity. Following previous work, we explored two 
types of abundance laws. As a first choice, we asumed that the 
molecular abundance is uniform across the core, and we used 
the Monte Carlo model to predict the intensity radial profile. 
Sometimes, this uniform abundance law predicted an emission 
profile that is much steeper than observed, and in these cases we 
modified the abundance law to include a central depletion hole 
of radius r,/ within which the abundance is negligible. 

As Fig.[5]shows, the N2H + radial profile of all the cores can 
be fitted using a uniform abundance model. Each core, however, 
requires a different abundance value. Cores A2 and C present the 
lowest N2H + abundances, with best-fit values that are four and 
two times lower than the value in core B, the most N2H + -rich 
core in the sample (see Table [4] for absolute abundance values). 
As mentioned above, the N2H + abundance in a core is expected 
to increase with time, so we can interpret the observed differ- 
ences in N2H + abundance as differences in the state of the chem- 
ical evolution of the cores. To better illustrate this evolution, we 
have re-ordered the cores in Fig. [5] so that the N2H + abundance 
decreases towards the right. If our understanding of the N2H + 
chemistry is correct, the sequence in the figure must represent 
a sequence of cores of decreasing age, in which core B is the 
oldest and cores A2 and C are the youngest. (Core D has been 
placed in the sequence by considering its N2H + (l-0) intensity 
towards the center, as no abundance estimate was possible with- 
out a mm-continuum map.) 

The analysis of the SO abundance provides further support 
to the N2H + -derived age sequence. As mentioned above, SO is 
expected to behave opposite to N2H + , and therefore show the 
largest abundances in the chemically youngest cores. Indeed, 
Fig. [5] and Table [4] show that the young cores A2 and C have 



large and uniform SO abundances, while the more evolved B 
and Al cores present SO depletion holes towards their center. 

Finally, all cores require a depletion hole in their C 18 abun- 
dance law to fit the observed emission profile. The youngest A2 
core requires the smallest hole, while the oldest B core requires 
the largest one (Table 0). Still, the correlation between evolu- 
tionary stage and C I8 depletion hole is weak, probably because 
determining the correct C 18 abundance law depends on mod- 
eling the poorly constrained envelope that surrounds each core. 
In any case, the results of the C 18 analysis are consistent with 
the interpretation that the L1517 cores are at different stages of 
evolution. 

Before finishing our analysis, it is worth mentioning that 
in the evolutionary sequence of Fig. [5] core D, for which no 
modeling could be carried out, appears as a relatively evolved 
core based both on its strong N2H + emission and on the flat- 
tened radial profiles of SO and C 18 0. This conclusion, however, 
disagrees with the recent analysis by iHirota et a l. (2009), who 
suggests that L1517D is an unusually young core from a chem- 
ical point of view. Resolving this contradiction requires making 
mm-continuum observations of this object to determine its true 
density structure and to allow a Monte Carlo modeling of its 
abundance profiles like that of the other cores. 

Although the cores in L1517 differ significantly in their 
chemical composition and therefore evolutionary state, they 
present very similar kinematic properties. Table[3]shows that the 
N2H + linwewidths are all close to 0.2 km s , which is typi- 
cal of the low-mass cores in Taurus and indicative o f subsonic 
nonthermal motions dMverslll983t ICaselli et ai]|2002l) . In addi- 
tion, the internal velocity gradients of the cores are on the or- 
der o f 1 km s" 1 pc" 1 , which is al so typical of low-mass cores 
(Goodma n et al.lll993t ICaselli et a l. 2002), and again suggests 
that the L 1 5 1 7 cores are representative of the population of low- 
mass cores in Taurus as a whole. This combination of differ- 
ent chemical composition and similar internal velocity structure 
suggests that the kinematic properties of the cores change little 
during their contraction from the more diffuse gas in the cloud. 
To further elucidate this issue, in the next sections we analyze 
the kinematics of the more extended gas in the cloud, and we 
compare it with that of the cores. 

5. Gaussian decomposition of the spectra 



In section I3TT1 we saw that each filament in the L1517 cloud 
is characterized by a well-defined range of LSR velocities 
(Table [TJ. The filaments, however, are not monolithic entities 
and seem to have an internal velocity structure characterized by 
low-level velocity gradients. As these gradients have potential 
information on the gas motions responsible for forming the em- 
bedded dense cores, we analyzed them in more detail by inspect- 
ing the individual C I8 O(l-0) spectra. From this inspection, we 
find that 92% of the spectra with peak intensity S/N> 3 present a 
single velocity component, while the rest of the spectra present 
evidence of two partly-overlapping velocity components. This 
preponderance of profiles with a well-defined number of com- 
ponents suggests that the velocity field at each cloud position 
can be characterized by fitting Gaussians to the spectra. 

To fit the more than 2,000 C 18 O(l-0) spectra observed to- 
wards L1517 we have used a semi-automatic procedure in the 
CLASS software and chosen to fit one or two components de- 
pending on the complexity of the line profile. The result from 
this fitting is illustrated in Fig[6] which also shows examples of 
the profiles. As can be seen, most of the cloud spectra have been 
fitted with single Gaussians (as shown in box number 4), and 
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Fig.5. Radial profiles of 1.2mm continuum, N2H + (l-0), SO(32-2i), and C 18 O(l-0) emission towards the dense cores in L1517. 
The cores were ordered by decreasing central emission of N2H + (l-0) to illustrate a possible evolutionary sequence (see text). The 
1.2mm data were averaged over 15"-wide intervals to increase legibility. In each panel, the solid red line represents the best fit- 
model, and when an additional dashed line is presented, the best fit required a central abundance hole and the dashed line represents 
the constant abundance model for comparison. No modeling of core D was attempted for lack of mm-continuum data. 



Table 4. Density and abundance profiles in the L1517 cores. 



(i) 





Density parameters 


Molecular abundances and 


depletion radii 


Core 


(cm- 3 ) 


r l/2 

(arcsec) 


a 


X(N 2 H + ) 


X(SO) 


r,,(SO) 
(cm) 


X(C 8 0) 


r rf (C 8 0) 
(cm) 


Al 
A2 

B 

C 


7.0x10* 
6.0x1 4 

2.2xl0 5 
4.7x10* 


60 
60 
35 
60 


2.5 
3.5 
3.5 

2.5 


1.0x10-'" 
0.4x10-'" 
1.5x10-'" 
0.7x10-'" 


l.OxlO- 9 
l.OxlO- 9 
0.4xl0- 9 
2.0x10-' 


1.3x10" 

0.0 
1.2x10" 

0.0 


1.5x10-' 
1.5xl0- 7 
1.5x10-' 

1.5x10-' 


1.5x10" 
1.1x10" 
1.9x10" 

1.7x10" 



Notes: (1) f All abundances are relative to H2. Core D was not modeled due to lack of mm-continuum data. 



the fits with two Gaussians are restricted to three well-defined 
regions identified with solid boxes. The largest of these three 
regions (box number 2) is centered near (Aa,A5) = (-200", 
+250") and seems to correspond to the superposition of fila- 
ments 3 and 4. This region appears in opti cal images a s a rel- 
ative enhancement in the obscuration, and iLee & Mversl ([1999) 
even classified it as an additional core (L1517B-2), but our ob- 
servations suggest that it is more likely a superposition effect. 
Another region with double spectra (box number 1) also seems 
to result from the superposition of two different components, 
this time filament 3 and an unlabeled and more diffuse compo- 
nent parallel to filament 4 that can be seen in the reddest map of 
Fig. |2 Finally, the third region with double spectra occurs near 
(Aa, A<5) = (-200", -200") and has a less clear origin. Its loca- 
tion would suggest an origin in the superposition of filaments 1 
and 2, but the velocity of the blue component does not match the 
velocity of any of the filaments. Observations with higher S/N 
towards this region are needed to determine the origin of this 
anomalous component. 



The Gaussian fitting procedure has also been applied to the 
SO and N2H + spectra. For SO, which is almost as extended spa- 
tially as C ls O (see Fig. 0J, we again find that we need to use 
two Gaussian components inside the regions discussed before, 
while single-peaked spectra are the norm in the rest of the cloud. 
For the N2H + (l-0) spectra, we fitted all hyperfine components 
simultaneously using the HFS method in CLASS and derived 
both the line centroid and the optical-depth corrected linewidth 
in this way. As shown in Fig.|4j the N2H + emission is much more 
compact than the C 18 and SO emission, and our inspection of 
the individual spectra found that all positions could be fitted with 
a single velocity component. 



The linewidths and velocity centroids derived with the 
Gaussian fits just described constitute the two input parameters 
for our analysis of the gas kinematics in L1517. In the follow- 
ing two sections, we analyze the behavior of each quantity sepa- 
rately. 



10 



A. Hacar and M. Tafalla: Fragmentation of velocity coherent filaments in L15 17 



a 

h o 

53 




Ac* (arcsec) 

Fig. 6. Graphical summary of the Gaussian fits to the C 18 O(l-0) 
emission in L1517. The gray scale shows the C 18 O(l-0) inte- 
grated intensity map of Fig. [T] and the superposed blue boxes 
enclose the regions where two Gaussians were fitted to the spec- 
tra. Insets with numbers 1, 2, and 3 illustrate double-Gaussian 
fits, while inset number 4 shows an example of the dominant 
single Gaussian fitted to most of the cloud spectra. 



6. Linewidth analysis 

6.1. Linewidth statistics: prevalence of subsonic motions 

The linewidth of a spectrum combines contributions from both 
thermal and nonthermal gas motions. In this section, we are in- 
terested in the nonthermal motions of the gas, which are poten- 
tially associated to gas turbulence or core formation motions. To 
isolate these motions, we have subtracted the thermal compo- 
nent from the measured linewidth following the standard prac- 
tice of assuming that the two contributions are indepen dent of 
each other so they add in quadrature (e.g.. iMversll 1983b . In this 
way, we estimate the nonthermal velocity dispersion for each 
species as 



cr NT - 



AV 2 kT 
81n2 ~ ~m' 



(4) 



where AV is the measured FWHM linewidth, k the Boltzmann's 
constant, T the gas kinetic temperature, and m the mass of the 
molecule under consideration. This velocity dispersion ct^t can 
be directly compared to the (isothermal) sound speed of the gas, 
c s , which has a value of 0. 19 km s -1 for ISM gas at 10 K. 

Figure Q shows the histograms of the velocity dispersion for 
the three molecular species in our survey (using all spectra with 
peak intensity S/N> 3). As can be seen, the overwhelming ma- 
jority of the spectra present a subsonic velocity dispersion and 
lie to the left of the cr^T lc s - 1 dashed line. The few positions 
that exceed the sonic limit form a low-level tail in the distribution 
that contains 2% of the spectra in C 18 and 3% of the spectra in 
SO. No N2FF spectrum has a supersonic linewidth. 

In addition to a very small fraction of supersonic points, the 
main feature of the histograms in Fig. [7] is the presence of a 
well-defined central peak. This peak indicates that the gas in 
the cloud is not only subsonic, but it also has a favored num- 
ber of nonthermal motions that approximately equal half the 
sound speed. More specifically, the histograms in the figure have 
the following mean and rms: (Tnt/c s = 0.54 + 0.19 for C 18 0, 
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Fig. 7. Histogram of the nonthermal velocity dispersion for C ls O 
(blue), SO (green), and N2H + (red) illustrating the predominance 
of subsonic values. 



cr NT /c s = 0.51 ± 0.23 for SO, and a NT /c s = 0.36 + 0.09 for 
N 2 H + . 

The values for the C 18 and SO nonthermal linewidth seem 
consistent with each other. The mean cr NT /c s of N2H + , on the 
other hand, is significantly smaller than that of C I8 and SO. 
Optical depth effects could potentially play a role in this differ- 
ence, since the hyperfine fit to the N2H + (l-0) spectrum auto- 
matically corrects for optical depth broadening, while no cor- 
rection is applied to the C I8 and SO data when fitting single 
Gaussians. To test whether optical depth broadening can explain 
the higher C I8 mean value, we have applied a simple correc- 
tion to the C 18 data by assuming an excitation temperature of 
10 K and estimating the optical depth from the peak intensity of 
the observed line profile. This correction is relatively small, as 
the optical depth never reaches unity, so even after its applica- 
tion, the distribution of cr NT still has a significantly higher mean 
value than for N2H + , and both a Kolmogorov-Smirnov and a 
Wilcoxon-Mann- Whitney test confirm that the two distributions 
must be different. 

The lower cr NT value of the N2I-F spectra could indicate that 
the dense core gas has an intrinsically lower velocity disper- 
sion than the less-dense ambient material and that core forma- 
tion has been accompanied by certain amount of dissipation of 
nonthermal motions. If this is the case, we can quantify the loss 
of nonthermal motions from the difference between the FWHM 
linewidths of C I8 and N 2 H + , which is 0.07 + 0.04 km s 1 . 

Even if small, the above estimate of the loss of non thermal 
motions during core formation is likely a significant overesti- 
mate. The nonthermal linewidth is not a local parameter, but the 
accumulated effect of gas motions along the line of sight. For 
this reason, some consideration should be given to line-of-sight 
effects when comparing tracers as different as C 18 and N2H + . 
As shown in section l4Tl the N2H + emission is significantly more 
concentrated than the emission from C I8 and SO because of its 
special chemistry and excitation requirements, so it samples a 
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shorter line of sight path than C 18 and SO. In addition, we see 
in the next section that part of the kinematics of the C 18 0/SO- 
emitting gas arises from large-scale motions in the filaments. As 
these motions likely accumulate along the line of sight, they are 
expected to contribute more to the C 18 and SO linewidth than 
toN 2 H + . 

To estimate the importance of this effect, we compare the 
sizes of the typical emitting regions for N2H + and C 18 0. Tables|2] 
and [4] show that the typical diameter of an L1517 core traced 
in N2H + is about 2 x 55", while the same parameter for a fila- 
ment traced in C 18 is about 2 x 95". In addition, Table|5]shows 
that the typical velocity gradient in the filament gas is about 
1.2 km s pc~'. Thus, if we assume a gradient of this size and 
multiply it by the difference in line of sight length between the 
two molecules, we predict a velocity difference between the trac- 
ers due to large-scale motions of about 0.06 km s . This differ- 
ence is so close to the observed linewidth difference that it seems 
very likely that a significant part of the nonthermal linewidth dif- 
ference between N2H + and C ls O (and SO) arises from their dif- 
ferent sampling of the cloud large-scale motions. Whether this 
effect can explain the full linewidth difference cannot be said, 
but it does show that the difference in local nonthermal motions 
(and thus dissipation during core formation) must be necessarily 
less than estimated before. 

A simple consequence of the subsonic nature of the nonther- 
mal gas motions in L1517 is the small role that they can play in 
supporting the cloud against gravity. The ratio between the non- 
thermal and thermal contributions to the gas pressure is given by 
Pnt/Pt = (o~nt/c s ) 2 , so even if we use the C 18 linewidth un- 
corrected for optical depth, we find that the nonthermal motions 
contribute to the gas pressure with only ~ 1/3 of the thermal 
value. This contribution is likely to represent an upper limit, be- 
cause the above pressure ratio assumes again that all nonthermal 
motions arise from a loc al "microscopic" turbu lent motions that 
increases cloud support ( Chandr asekharll 1 95 lb and ignores the 
contribution from large-scale flows. 

6.2. Spatial distribution and consequences for turbulent 
models of core formation 

To further characterize the nonthermal motions, we now study 
their spatial distribution. Figure [8]presents a plot of cr NT /c s for 
C 18 and N2H + as a function of position along each of the fila- 
ments identified in section [3~T1 In agreement with the histogram 
analysis, all four filaments are dominated by subsonic gas over 
their entire length, and only a few positions exceed the sonic 
limit in C 18 0. These few supersonic positions (some of them 
coincident with filament overlaps) appear almost randomly scat- 
tered over the filaments length, and show no correlation with 
either the position of the cores or the beginning or end of the 
filaments. 

In addition to a lack of supersonic points, the plots of non- 
thermal linewidth in Fig. [8] present very little spatial structure. 
Filaments 1, 3, and 4 have almost flat distributions of C 18 
ctnt/c and well-defined mean values close to 0.5 (Table [5]). 
Filament number 2 also has a mean (T^tIc close to 0.5, but 
presents changes in the linewidth along its length. These changes 
are of unclear origin, and seem to result from the presence of two 
localized regions of enhanced (but mostly subsonic) linewidth. 
One of these regions corresponds to the vicinity of core B, and 
the other occurs in the vicinity of core A 1 . Concerning the N2H + 
linewidths, Fig.[8]again shows their tendency to be smaller than 
those of C 18 0. This trend is present in all the cores, although it 
is somewhat reduced in core D of filament 3. Apart from this 
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Fig. 8. Distribution of nonthermal velocity dispersion as a func- 
tion of position along each of the L1517 filaments. Blue squares 
represent C O(l-0) and red triangles and N2H + (l-0) values. 
The spatial scale for filament 3 has been shrunk by a factor of 
1 .5 in order to fit it into the reduced box size used for the other 
filaments. The mean formal error in the <tmt fc s determination is 
0. 1 for both N2H + and C I8 0, which is approximately the marker 
size. 

Table 5. Kinematic properties of the filaments (1) 



Filament 


(cr NT )/c s 


(Vlsr) 
(kms- 1 ) 


<|WiM|> 

(km s~'pc~') 


1 

2 
3 
4 


0.57+0.15 
0.63+0.17 
0.53+0.16 
0.41+0.16 


5.52+0.07 
5.79+0.07 
5.89+0.05 
6.14+0.04 


1.0+0.5 
1.4+0.7 
1.3+1.0 
0.9+0.5 



Notes: (1) From C 18 O(l-0) emission. 



already discussed effect, Fig.[8]shows that the N2H + linewidths 
remain approximately constant inside each of the cores. 

The flat distribution of the velocity dispersion in the fil- 
aments, and in particular in the vicinity of the dense cores, 
poses a significant constraint to models of core formation by 
supersonic turbulence. In these models, the cores are shock- 
compressed structures formed at the stagnation point between 
two convergent flows, and as a result, they are expected to 
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Fig. 9. Comparison between the C 18 nonthermal velocity dis- 
persion around each L1517 core and the predictions from a 
model of turbulent core formation. Top: schematic view il- 
lustrating how the velocity dispersion has been averaged inside 
four concentric shells around each core. Bottom: comparison 
between the observed C 18 data (r ed solid triang l es) an d the 
predictions from the LSD model by Klessen et al. (2005). The 
black solid squares represent the three different projections of a 
fully subsonic core, while the open inverted triangles represent 
the three projections of a core that has transonic motions in one 
of them. 



be surrounded (even c onfined) by turbulent layers of cloud 
gas (e.g., Padoan et al. 2001; Ballesteros-Paredes et al. 2003; 
iKlessen et al.ll2005l) . In L1517, the observed C 18 emission sur- 
rounds the location of the dense cores and arises from gas that 
is about one order of magnitude less dense than the core cen- 
ters (section l4Tb . so we can naturally expect that the C 18 emis- 
sion arises from the gas layers that physically surround the dense 
cores. (C 18 does not sample gas inside the cores due to freeze 
out, see section [4]) If the cores have formed by the convergence 
of gas flows, part of the C I8 emission must originate in gas 
in the flows and must retain a signature from the core-forming 
convergent motions. To test whether such a signature is present 
in our data, we calculated for each core the average C ls O ve- 
locity dispersion inside rings around the core center having radii 
equal to 1,2, 3, and 4 core radii (core radius is n/2 in Table @), 
This procedure is illustrated in the upper panel of Fig. [9] and has 
been stopped at four core radii to ensure that we sample equally 
gas that surrounds the cores in all directions, and not only in the 
direction of the large-scale filaments (which dominates the emis- 
sion at larger radii). The results from this average for each of the 
five cores in L1517 are indicated in the bottom panel of Fig. [9] 
Not surprisingly given the general trend seen in the filaments, 
all cores present a distribution of velocity dispersion with radius 
that has no systematic change with distance from the core center. 
To compare the distribution of velocity dispersions around 
the cores of LI 5 17 with the predictions from turbulent mod- 
els, we chose the numerical simulations of gravo-turbulent core 



formation bv lKlessen et al.l (12005). These authors attempt to re- 
produce conditions similar to those of a low-mass star-forming 
region like L1517 and present a very complete view of the in- 
ternal kinematics of the simulated cores. From the two families 
of simulations that these authors present, we selected the large- 
scale driving (LSD) case because it produces the most quies- 
cent co res, and is the refore more likely to fit our observations of 
L15 17. IKlessen et a l. (2005) present in their Fig. 1 radial pro- 
files of the velocity dispersion for three different views of two 
selected cores (likely chosen for their similarity to observations), 
and we used these profiles to derive nonthermal velocity compo- 
nents by subtracting in quadrature the contribution from thermal 
motions (as done with the C 18 data). These data have been av- 
eraged as a function of core radius by converting the normalized 
column densities given by the authors into a distance from the 
core center by assuming that the the density of the cores follows 
the profile of a Bonnor-Ebert sphere (as claimed by the authors). 

The pred i cted n onthermal velocity dispersions from the 
IKlessen et alj (2005) model are represented in Fig. [9] As ex- 
pected, these predicted dispersions show a significant increase 
with radius due to the presence of converging flows around the 
cores. The predicted increase is approximately a factor of 2 be- 
tween 1 and 4 core radii, and as a result, the outermost points 
in the graph are expected to exceed the sound speed limit in all 
core models. Even in the most quiescent, subsonic core model, 
which is arguably not representative of the sample because sub- 
sonic cores are less than 25% of the total in the simulation, the 
match between model and data is only acceptable for the inner- 
most two radii. For larger radii, model and data diverge with in- 
creasing distance from the core center, and end up differing by a 
factor of 2 in the outermost layers sampled by our observations. 
Even larger disagreement occurs for the case of transonic cores, 
which constitute ~ 50% of the total cores in the simulation and 
are therefore more representative of the model results. As can be 
seen, two out of the three sets of points exceed the sound speed 
at all radii, and therefore fail to fit the observations at all radii in 
all the cores. 

Although limited, the above comparison illustrates the basic 
disagreement between our observations of the L1517 cloud and 
the predictions from models of core formation by convergence 
of supersonic flows: the C 18 spectra in L1517 sample the core 
outer layers far enough to trace any core-forming motions, so 
if these motions were supersonic, they should have left a clear 
signature in the C 18 spectra. The absence of such a signature 
rules out the presence supersonic motions around the cores and 
sets a limit to any core-forming gas flow that is well inside the 
subsonic regime. 

6.3. Velocity coherent filaments 

In their st udy of the linewidth of dense-gas tracers in cores, 
Barranco & Goodman ( 1998) and Goodman et al. ( 1998) found 
that the nonthermal component remains almost constant over 
the core interior and that this behavior represents a departure 
from the well-known line width-size rel ation commonly associ- 
ated to turbulent motions (Larson 1981). These authors refer to 
this property of the core gas as velocity "coherence," and suggest 
that the scale of coherence, which they found to approximately 
coincide with the core size, may indicate that core formation is 
related to the process of turbulence dissipation. 

Our analysis of the C 18 linewidth illustrated in Figs.[8]and 
[9] shows that the region with constant nonthermal linewidth can 
be followed in significantly larger scales than the cores and that 
it extends over distances as long as the filaments themselves (~ 
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0.5 pc). This coherence of the velocity field in the filaments is 
not limited to the nonthermal linewidth, and it can be seen in 
the behavior of the velocity centroid, to be studied in detail in 
the next section. In the LI 5 17 clou d, therefore, the filaments are 
veloc ity coherent in the se nse of iBarranco & Goodman (1998) 
and iGoodman et alj (1 19981) . and this extends the scale size of 
coherence in the gas by a factor of approximately 5. 

Finding velocity coherence on scales as large as the filaments 
separates the scale of coherence from that of the cores, and sug- 
gests that core formation is not likely the direct result of turbu- 
lence dissipation. As seen in section |6T1 some turbulence dis- 
sipation may have occurred during core formation, but by being 
significantly subsonic, it is unlikely to have affected the gas pres- 
sure balance and triggered core formation. Velocity coherence in 
cores, therefore, does not seem a defining property of the con- 
densations, but a condition inherited from the larger spatial scale 
of the filaments. 

Our finding of velocity coherence on scales larg er than a core 
would seem to contradict the recent results from iPineda et alJ 
d2010l) . who claim to have detected of the transition between the 
coherent and turbulent regimes in a core with observations of 
the B5 region in Perseus. We note, however, that although these 
authors associate the velocity coherent region with a dense core, 
one can see in their Fig. 3 that the velocity coherent region in 
B5 is elongated and has a length of about 0.5 pc, similar to our 
L1517 filam ents. The B5 regio n, in addition, contains multiple 
dense cores (lEnoch et alj|2006i also our own unpublished data), 
and this suggests again that it is more similar to the L1517 fila- 
ments than to an isolated dense core. Clearly more observations 
of different core environments are needed to understand the size 
and relation between the different velocity coherent regions. If 
the L1517 cloud is a representative region of dense core forma- 
tion, we can predict that filament-wide velocity coherence will 
be a common phenomenon. 

7. Line centroid analysis 

The second output from our Gaussian fit to the spectra is the 
velocity centroid. Figure \TU\ compares the SO and N2H + cen- 
troids with those of C 18 towards all positions where the fits 
were considered to be significant (S/N > 3). As can be seen, 
the velocities of the different tracers agree with each other in- 
dependently of any variations in the bulk velocity of the cloud. 
The average difference in velocity between C ls O and either SO 
or N2FF is 0.03+0.03 km s _1 , which means that the velocities 
of the different tracers differ on average by less than one fifth 
of the sound speed. This good match between tracers rules out 
any significant motions between the different density regimes of 
the gas and, in particular, it rules out any systematic drift be- 
tween the dense cores (traced by N2H + ) and the surrounding gas 
(traced by C 18 0). S uch a quiescent sta te of the gas is not pe- 
culiar to L1517, and lWalsh etafl (12004 and lKirket ai1(l2007ah 
have found a similar lack of velocity shifts between N2FF and 
C 18 in a number of low-mass star-forming regions. 

7. 1. Continuity of the velocity field and large-scale 
oscillations 

The lack of velocity shifts between tracers does not mean that 
the gas in L1517 is static. Figure [TUl shows how the LSR veloc- 
ity of the cloud material spans almost 1 km s -1 , and this velocity 
spread arises from a combination of the different velocities of 
each of the filaments and the presence of internal velocity gra- 
dients inside the filaments. To more fully study these gradients, 
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Fig. 10. Comparison between the Vlsr of C 18 and SO (blue 
squares), and C I8 and N2FF (red triangles) in L1517 illustrat- 
ing the good agreement between all tracers. The central dashed 
line indicates the locus of equal velocities, and the surrounding 
dotted lines have been displaced by the sound speed (0. 19 km s" 1 
for ISM gas at 10 K). The mean formal error in the Vlsr deter- 
mination is 0.02 km s" 1 for SO and 0.01 km s _I for both C ls O 
and N2FF, which is comparable or smaller than the marker size. 



we present plots of the velocity centroids of C I8 and N2H + as 
a function of position along each filament in Fig.fTTI 

As can be seen in Fig.QT| the velocity field of each filament 
presents a remarkably low level of spatial change. The dispersion 
of the C I8 velocity centroids does not exceed 0.07 km s _I in 
any filament, even when considering a few regions of enhanced 
dispersion that seem to coincide with locations of multiple- 
peaked spectra and to overlap between different filaments. In ad- 
dition, the end-to-end velocity change in each filament is about 
0.2 km s or less. This quiescent state of the gas on scales as 
large as 0.5 pc once more justifies the interpretation of the fila- 
ments as velocity-coherent structures. 

Although the changes in velocity along the filaments are 
small, the low dispersion of the centroid data allows the exis- 
tence of large-scale velocity patterns to be discerned in some 
of the filaments. Figure QT] shows how the C I8 centroids os- 
cillate quasi-periodically, especially in filaments 1 and 2, over 
the length of the filaments with typical wavelengths of approxi- 
mately 0. 1-0.2 pc. As also shown in the figure, the more spatially 
localized N2FF centroids follow the C I8 oscillatory pattern at 
those positions where both tracers can be observed (i.e., towards 
the dense cores), indicating that both species trace the same gen- 
eral velocity pattern. This good match between the C ls O and 
N2I-F data is remarkable because the two species do not coexist 
spatially due to their antic orrelated chemistries (section @), and 
they therefore trace significantly different regimes of gas density. 
Thus, the continuity of the C 18 and N2FF velocity gradients 
must reflect a continuity between the velocity field of the fila- 
ment gas and the internal velocity gradients of the dense cores 
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Fig. 11. Velocity centroid as a function of position along each of 
the L1517 filaments. Blue squares represent C 18 O(l-0) and red 
triangles N2H + (l-0) values. The spatial scale for filament 3 has 
been shrunk by a factor of 1.5 in order to fit it into the reduced 
box size used for the other filaments. The mean formal error in 
the velocity centroid determination is 0.02 km s for C I8 and 
0.01 km s _1 for N2H + , which is always smaller than the marker 
size. 



embedded in it. This implies that the internal velocity gradients 
in the cores are not intrinsic core properties (like isolated ro- 
tation), but result from the large-scale motions of surrounding 
filament. Core velocity gradients should therefore be interpreted 
in terms of the velocity gradients of the lower density gas. 

The observed continuity between the velocity field of the fil- 
aments and the cores implies that core formation in L1517 has 
not decoupled the dense kinematically gas from its surrounding 
environment. This behavior again seems to contradict the expec- 
tation from the turbulent models of core formation, because in 
these models, the cores are formed by the direct shock of streams 
of lower density gas, so a discontinuity in the velocity field is ex- 
pected. The observations of L 15 17 are better understood if the 
filaments are coherent in velocity as a whole, and the cores have 
formed from the contraction of this velocity-coherent gas. The 
transition from cloud to core conditions must have therefore in- 
volved little or no dissipation of kinetic energy, and in particular 
the absence of supersonic shocks. The data therefore indicates 
that core formation in L1517 has been an almost quasi-static 
process. 



7.2. A simple model of the velocity oscillations 

The velocity oscillations of filaments 1 and 2, dominate the 
large-scale kinematics of the gas in these objects, and suggest 
some type of ordered pattern in their three-dimensional veloc- 
ity field. Our limited observations cannot fully constrain the 3D 
properties of this velocity field, but the large-scale nature of the 
motions and the lack of systematic velocity gradients perpen- 
dicular to the filament axis suggest that the underlying velocity 
field must have an important component directed along the axis 
of the filaments. If this is correct, it is intriguing to investigate 
whether the velocity oscillations of filaments 1 and 2 are related 
to the presence of embedded cores, and in particular, whether 
they could represent part of the motions responsible for core for- 
mation. 

Core formation in a filamentary gas cloud is often in- 
vestigated using the idealized geometry of an infinitely-long, 
axially-symmetric gas cylinder. The stability of such a 
configuration has been studied in detail u s ing both semi- 



analytical tec h niques ([St odolkie wiczl 1963b lOstriker l 
Larsonj_ 1985b iNagasawal 1987b Inuts uka & Mivama 



1964 



1992; 



Nakamura'et al.1 ll993rTGehmanet alj 11996b iFiege & Pudritz 



2000aj) and numerical simula tions (Bastien 1983; Bastien et al. 



1991b iNelson & Papaloizoul 1 1993b INakamura et all 1 1995b 



iFiege & Pudritzl l2000bl) . From this body of work we know 
that, depending on the initial equilibrium state of the gas, the 
amount of support provided by pressure and magnetic fields, 
and the geometry of both the magnetic field and the perturbation 
applied to the system, the gas cylinder can follow a number of 
evolutionary paths, including expansion, collapse to a spindle, 
or fragmentation into multiple clumps. Among these outcomes, 
cloud fragmentation is the most interesting here because of its 
possible relation to core formation and the velocity oscillations 
in filaments 1 and 2. 

In the simplest case of an isothermal cylinder in equilibrium, 
fragmentation occurs through a Jeans-type of instability in which 
a density perturbation with a large-enough wavelength accumu- 
lates enough mass in each over-dense region to make it gravita- 
tionally unstable (e.g., Stodolkiewicz 1963; Larson 1985). The 
most unstable perturbation is an axisymmetric mode that breaks 
the initially continuous filament into a chain o f dense conden- 
sation s equally spaced along the cylinder axis ( INakamura et al.l 
1993). Numerical simulations show that this process involves 
redistribution of the gas in the filament via motions that have 
a dominant velocity component parallel to the filament axis, at 
least during the first stages o f evolution (e.g., Nakamur a et al.l 
Il993b IFiege & P udritz 2000b). The similarity of these motions 
with those inferred for filaments 1 and 2 in L1517 is the main 
motivation for attempting a simple kinematic model. 

Modeling in detail the fragmentation of the L1517 filaments 
exceeds the scope of this paper, and would require informa- 
tion on so far unknown cloud properties like the strength and 
geometry of the magnetic field. For this reason, we limit our- 
selves to testing whether the velocity oscillations of Fig. QT| are 
consistent with large-scale core-forming motions, and we do so 
with a highly idealized infinite cylindrical model. Our first ap- 
proximation is to only study motions along the central axis of 
the cylinder, which allows us to convert the problem to one di- 
mension and which is motivated by the dominance of longitu- 
dinal motions in sim ulations. Following standard practice (e.g., 
Binn ev & Tremainel 1987. their section 5. 1), we use perturbation 
analysis and assume that the gas starts from a state of hydrostatic 
equilibriu m in which the density is constant along the filament 
axis (e.g., Ostriker 1964). To this equilibrium state we added a 



A. Hacar and M. Tafalla: Fragmentation of velocity coherent filaments in L15 17 



15 




1 i 1 r 




CD 

E 

CO 



100 200 300 400 500 600 
L (arcsec) 



C) 



E 
> « 






> 

CO 




100 200 300 400 500 600 
L (arcsec) 

Fig. 12. Kinematic modeling of the velocity oscillations in filaments 1 and 2. a: Schematic view showing how a core-forming 
velocity field along the filament axis causes a A/4 shift between the sinusoidal perturbations of density and velocity, b and c: 
Comparison of the density and velocity perturbations in filaments 1 and 2 with the expectation from the simple kinematic model. 
For each filament, the upper panel (black squares) shows the increment of mm continuum flux over the mean (a proxy for the density 
perturbation), and the lower panel (blue squares) shows the variation in the C 18 centroid velocity over the filament mean (after 
subtraction of a linear gradient). The green solid line in each upper panel is a sinusoidal fit to the density perturbation, from which 
a wavelength and phase are determined. The green solid line in each bottom panel is the result of shifting the density sinusoid by 
A/4 (plus scaling it arbitrarily), and represents the expected pattern for a core-forming velocity field. The mean formal error in the 
SVlsr determination is 0.01 km s for both filaments, which is on the order of the marker size. 



small longitudinal perturbation in both density and velocity, and 
work only to first order in the perturbation quantities (linear anal- 
ysis). Under these conditions, the equation of continuity can be 
written as 



dpi du\ 







(5) 



where p is the density, u the velocity, z the spatial coordinate 
along the filament axis, and the and 1 subscripts refer to the 
unperturbed solution and the first-order perturbation. If we now 
assume Fourier-component perturbations for both density and 
velocity, we can write 



Pi(z,t) = C\ e 



i(kz-u)t) 



and ui(z, t) - Ci e' 



i(kz-u)t) 



(6) 



where C\ and C2 are constants due to the lack of z-dependence of 
the unperturbed solution. Substituting these expressions into the 



equation of continuity, we can obtain a simple relation between 
the two constants. The case of interest here is that of an unstable, 
core-forming mode (i.e., opposed to a stable sonic wave), and 
this implies that to 2 must be negative. Defining a real number y 
so lo - iy, we can write 



r .Cxy 
C 2 = i — -. 
Po k 



(7) 



As po, y, and k are real numbers, the above equation implies 
that the Fourier component of the velocity perturbation must 
be shifted in phase from the density component by n/2, which 
is one quarter of th e wavelength of the perturbation (see also 
Geh man et al.ll 19961) . (In a stable (sonic) perturbation, a> is a real 
quantity, so the density and velocity Fourier components are in 
phase.) 

The origin of the A/4 shift between the velocity and density 
perturbations can be easily understood with the simplified car- 
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toon of Fig. Q~2]a. As can be seen, for the gas motions to be 
core-forming, they have to converge towards the core centers, 
and this means that a density peak must correspond to a position 
of vanishing velocity. Assuming that both density and velocity 
perturbations are sinusoidal, this requires that there is a A/A shift 
between the two. 

The expected A/A shift between the density and the velocity 
patterns in an unstable perturbation provides a simple criterion 
to test whether the velocity field in the L1517 filaments is con- 
sistent with core-forming motions. To this end, we have fitted a 
sinusoid to a cut of the mm-continuum emission along the fila- 
ment axis. This quantity should be proportional to the filament 
density profile, with the caveat that the bolometer observation 
filters the extended emission due to chopping, and from this fit 
we have determined the wavelength and phase of the density per- 
turbation in filaments 1 and 2. The top panels of Fig.Q~2]b and c 
show that both filaments can be reasonably fitted with sinusoidal 
density profiles that have wavelengths of approximately 340 and 
270 arcsec, respectively. 

To now test whether the observed velocity field is consistent 
with the formation of these density profiles, we shift the density 
sinusoid by one quarter of the wavelength, as required by the 
perturbation analysis, and compare the result with the velocity 
centroid data towards the filament axis. For this step, we have 
one degree of freedom in the choice of the sign of the shift, as 
there are two possible inclinations of the filament with respect to 
the plane of the sky (e.g., the eastern part of filament 1 could be 
inclined toward us or away from us), and this adds a sign ambi- 
guity to the gas radial velocity. When this inclination is chosen 
(together with a scaling of velocity amplitude), we derive the 
model predictions shown in the bottom panels of Fig.[T2]b and c 
superposed to the observed velocity centroids (to which we have 
subtracted global linear gradients of 0.5 km s pc~' in filament 
1 and 0.1 km s -1 pc~' in filament 2). 

As can be seen, the velocity field of filament 1 agrees rea- 
sonably well with the prediction from the shifted sinusoid both 
in wavelength and phase, suggesting that the observed velocity 
oscillation is consistent with core-forming motions. The velocity 
towards core C does indeed coincide with a zero value in the ve- 
locity sinusoid, while the behavior of the velocity towards core 
A2 is not as clear due to the higher scatter and the few velocity 
points. Such a reasonable fit is encouraging for an interpretation 
of the velocity field in terms of fragmentation, as this filament is 
also the best behaved according to both the maps (Fig. [2]) and the 
radial profiles (Fig. |3). It also harbors the chemically youngest 
cores of the cloud (Fig. [5}, suggesting that it is at the earliest 
phases of fragmentation. 

In contrast to filament 1, filament 2 does not fit the shifted 
sinusoid velocity pattern as well. On the one hand, the veloc- 
ity field of this filament does not follow a sinusoid pattern very 
closely, and on the other, a region of high scatter between cores 
B and Al leaves the velocity poorly defined, especially in the 
vicinity of core Al. Still, as the figure shows, the position of 
core B lies in a region consistent to have zero velocity, which is 
the expected pattern for core-forming motions. 

So far we have concentrated our analysis on the phase of the 
velocity oscillation, but the amplitude can also be measured and 
contains information on the speed of the possible core-forming 
motions. As illustrated by the cartoon in Fig. |T2j the observed 
amplitude represents only the line-of-sight component of the true 
velocity amplitude, and it therefore needs to be corrected for the 
angle that the filament makes with the line of sight. Lacking a 
better estimate, we assume that each filament makes a 45 de- 
gree angle with the line of sight, and we use this value to de- 



rive representative velocities of the possible core-forming mo- 
tions. As can be seen from Fig.[T2l both filaments have velocity 
oscillations with an amplitude of approximately 0.04 km s -1 , 
which would correspond to a projection-corrected value of ap- 
proximately 0.06 km s . I nterestingly, such subsonic speeds are 
similar to those inferred bv lLee et al.l (11999) for the inward mo- 
tions towards a large number of low-mass dense cores, which 
were obtained using a very different technique (analysis of self- 
absorbed profil es). Also in agreem ent with previous core life- 
time estimates dLee & Mverslll999l) . the time scale of collapse 
inferred from the y value in Eq.|6]is on the order of 0.5 Myr. 

Even if the velocity oscillations in L1517 do not arise from 
contraction motions, their amplitude constrains core-formation 
models, as any underlying contraction motion has to be slower 
than about 0.05 km s -1 to remain undetected or confused with 
the observed oscillation. Thus, no matter how we interpret the 
observed velocity pattern, the conclusion that core formation in 
the L1517 cloud is strongly subsonic seems inescapable. 

8. How did the L1517 cores form? 

The picture that emerges from our analysis of L15 17 is that core 
formation in this cloud has involved only subsonic motions and 
that it has started from conditions that were very close to hydro- 
static equilibrium in an elongated configuration. These charac- 
teristics, together with the close fit of the filament density pro- 
files by the model of an isothermal pressure-supported cylinder, 
suggest that some form of cylindrical gravitational fragmenta- 
tion has played a role in the formation of the LI 5 17 cores. 

As discussed in section I7T21 the simplest case of cylindrical 
fragmentation is that of an isothermal, infinitely-long, pressure- 
supported cylinder. Perturbation analysis shows that the gravita- 
tional fragmentation of such a system has a critical wavelength 
Af — 3. 9 AH, where H is the filament width given by equation [2] 
(Stodolkie wiczl [19631) . Core formation via this simplest gravi- 
tational fragmentation model therefore requires that the dense 
cores are physically spaced by a distance grater than Af. 

To test whether the separation between cores in L1517 is 
consistent with the simplest gravitational fragmentation sce- 
nario, we have estimated the critical wavelength of the filaments 
using the H values obtained from modeling their radial profiles 
(Table |2). These values imply that Af is about 1.2 x 10 18 cm, 
or 570" for our assumed distance to Taurus (144 pc). As can be 
seen from Fig. Q~2] the observed distance between the cores in 
filaments 1 and 2 is « 300", which is significantly less than the 
minimum value expected from the theory. 

Assuming that the filaments are inclined with respect to the 
plane of the sky can change the above estimate, but it does not 
bring significant improvement. On the one hand, correcting for 
projection will increase the true distance between the cores, but 
it will also increase Af because the projection correction af- 
fects the column density estimate from which Af is derived (e.g., 
Arzoumanian et al. 2011). As a result, the observed separation 
between cores in the L1517 filaments seems to be about 1.5- 
2 smaller than predicted by the simplest model of gravitational 
fragmentation. 

A number of factors can potentially explain the smaller- 
than-predicted separation between cores in L1517. First of all, 
our Af value could have been overestimated. This is possible 
because Af was derived from the density analysis of the fil- 
aments, which depends on assuming a standard C I8 abun- 
dance. If this molecule has suffered from large-scale freeze out 
in the filaments, the true gas volume density would be higher, 
and Af correspondingly smaller. A second possibility is that 
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edge effects, due to the finite size of the filaments, are impor- 
tant in the outcome of the fragmentation, as suggested by the 
simulations of Bastien ( 1983) and Nelson & Papaloizou (1993). 
Unfortunately, these effects depend very strongly on how sharp 
the filament edges are ( Nelson & Papaloizou 1993), and our data 
are not detailed enough to constrain this parameter. Finally, mag- 
netic fields are known to affect the fragmentation of a cylinder 
(Nakamura et al.ll993l [T995; Geh man etail 1996; Hanawa et al. 
1993; Fiege & Pudritz 2000b). Fiege & Pudritz (2000b), in par- 
ticular, show that both toroidal and poloidal magnetic fields de- 
crease the critical length of fragmentation, and can do so by a 
significant factor for sufficiently strong fields. (They also de- 
crease the instability growth rate.) Whether this could explain 
the observations of L1517 is u nclear, as l i ttle is known about 
the magnetic field in the cloud. iKirk et all (2006) find that the 
field towards core B is too weak to provide support, but these 
measurements are limited to a small region of the cloud already 
known to have undergone core formation. A detailed character- 
ization of the large-scale magnetic field in the L1517 cloud is 
clearly needed. 

Even if the L1517 cloud is less symmetric than assumed 
by our model, and its fragmentation history more complex than 
predicted by the simple theory, the observations presented here 
clearly constrain the process of core formation to have involved 
an elongated geometry and mostly subsonic motions. Core for- 
mation in this environment therefore appears to be a two-step 
process in which the elongated configuration close to equilib- 
rium is formed first and subsonic fragmentation into cores oc- 
curs later. This two-step scenario of core formation has already 
been proposed by several authors based on different considera- 
tions, and it therefore seems to apply to more regions than just 
the one studied he re (|Schneider & Elmegr eenlll979l; iHarfmannl 
l2002tlMvers1l2009l) . What our L1517 data show now is that the 
filamentary gas prior to fragmentation has a subsonic level of 
turbulence and a coherent velocity pattern, two properties that 
were thought before to only apply to dense cores. For such a 
quiescent state of the gas to occur, turbulent motions must have 
dissipated prior to (or rapidly during) filament formation, and 
therefore must play a limited role in the formation of the indi- 
vidual cores. If such a scenario is representative of core forma- 
tion in other dark clouds, subsonic, velocity-coherent filaments 
must be a common feature, and embedded cores must systemati- 
cally present kinematic patterns continuously connected to those 
of the large-scale filaments. Detailed study of the connection be- 
tween dense cores and their surrounding cloud in a larger sample 
of regions should be able to test these two predictions. 

9. Summary 

We observed the L1517 dark cloud in C I8 O(l-0), N 2 H + (l-0), 
SO(Jn=32-2i) with the FCRAO telescope, and in the 1.2mm 
dust continuum with the IRAM 30m telescope. From the analy- 
sis of these data, we came to the following main conclusions. 

1. The gas in the L1517 cloud is structured in four filaments 
with typical sizes of about 0.5 pc. 

2. The radial profile of C 18 O(l-0) emission in each filament 
consists of a central flattened region and a power-law tail. An 
isothermal cylinder model provides an approximate fit to the 
emission, although better agreement is obtained with a softened 
power law profile. 

3. Five starless cores are embedded in the filaments. Their 
chemical composition indicates that they are at different evolu- 
tionary stages, although their kinematic properties are very sim- 
ilar. 



4. The velocity field of the gas in the L1517 cloud has been 
characterized by fitting Gaussians to the observed spectra. Most 
positions require one Gaussian, although a few regions require 
two, very likely due to the overlap between different compo- 
nents. 

5. The filaments are extremely quiescent. Their nonthermal 
linewidth is subsonic and changes very little over the length of 
the filaments. The velocity centroids also change subsonically 
over the filament length. These characteristics indicate that the 
gas in the filaments is velocity coherent on scales of about 0.5 pc. 

6. Although quiescent, the filaments have large-scale pat- 
terns of (subsonic) velocity. The gas in the dense cores follows 
the velocity pattern of the less dense gas closely, indicating that 
core formation has not decoupled the dense gas kinematically 
from its surrounding filament material. 

7. In two filaments, the large-scale velocity patterns seem to 
consist of oscillations. A simple kinematic model shows that, at 
least in one filament, the oscillatory pattern is consistent with 
core-forming motions along the axis of the filament. 

8. Core formation in L1517 seems to have occurred in two 
steps. First, the subsonic, velocity-coherent filaments have con- 
densed out of the more turbulent ambient cloud. Then, the cores 
have fragmented quasi-statically and inherited the kinematics of 
their parental filament. Turbulence dissipation has therefore oc- 
curred on scales of 0.5 pc or larger, and seems to have played 
little role in the formation of the individual cores. 
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